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Abstract 

A phonon-roton dispersion relation is proposed for the elementary 
excitations of a quantum biexcitonic liquid in semiconductors. The pro- 
posed dispersion relation is used as a starting point for a calculation 
of the photoluminescence spectrum of the liquid and an analysis of its 
behavior under variation of the temperature and density of the biexci- 
tonic liquid. The parameters of the dispersion curve of elementary ex- 
citations of the quantum biexcitonic liquid are evaluated by fitting the 
calculated photoluminescence spectrum to the experimental spectrum of 
the biexcitonic liquid of semiconducting /3-ZnP2 crystals. Experimental 
studies of how the photoluminescence spectrum of a biexcitonic liquid 
in /3-ZnP2 depends on the temperature and the intensity of the laser 
excitation confirm the initial theoretical model. The dependence of the 
temperature of the crystals on the excitation intensity is measured, and 
for some of the samples an anomalous dependence is found: the tem- 
perature of the crystal decreases as the excitation intensity increases. 
This effect is probably a consequence of the giant thermal conductivity 
of the superfluid biexcitonic liquid in /3-ZnP2 crystals. 

Keywords: biexcitons, quantum biexcitonic liquid, phonon- 
roton 

dispersion relation, superfluidity 
PACS: 71.35.Lk, 78.55-m 



Deceased 



1 



Introduction 

For a long time the phenomenon of superfluidity was observed 
exclusively in liquid '^He at temperatures below 2.17 K. It is well 
known that the dispersion relation for elementary excitations in 
superfluid He II is characterized by a so-called roton gap. A dis- 
persion relation of this sort was first proposed by Landau [jlj and 
later investigated theoretically by Feynman,[]2[ and it was subse- 
quently confirmed in neutron scattering experiments.]^ Among 
the Bose liquids, ^He is unique in that it remains liquid at atmo- 
spheric pressure down to the lowest temperatures. Other Bose 
hquids crystalhze at temperatures above those at which superflu- 
idity can arise. Crystallization sets in because the energy of the 
intermolecular interaction is greater than the vibrational energy 
at temperatures above the temperature of the onset of superflu- 
idity. Since the vibrational energy goes as /icjvibr oc m^^/^ (m is 
the mass of the molecule), superfluidity is more likely to arise in 
liquid molecular hydrogen than in other liquids. However, even in 
the case of liquid H2 the energy of the intermolecular interaction 
is too large, and crystallization sets in before superfluidity as the 
temperature is lowered. 

The analogy between Wannier-Mott excitons in semiconduc- 
tors and atomic hydrogen is well known. This analogy extends 
to excitonic molecules (EMs or biexcitons) and to the hydrogen 
molecule. The problem of Bose condensation and superfluid- 
ity in a system of excitons has been very popular in semicon- 
ductor physics for the past 20-25 years (see, e.g., Refs.[^-[^). 
The possibility of superfluidity of excitons in a strong magnetic 
field was pointed out in Ref.[]5j. An unusual ballistic solitonlike 
regime of exciton motion in CU2O crystals was observed in Ref.[g]. 
This effect was interpreted as evidence of superfluidity of Bose- 
condensed excitons. 

Guided by the analogy between excitonic molecules and hy- 
drogen molecules, one can analyze the possibility for superfluidity 
to arise in such a biexcitonic (molecular) liquid. The existence 
of a biexcitonic liquid in semiconductors was first considered in 
Ref.[|7j. As we have said, superfluidity is not observed in liquid 
hydrogen because of the earlier onset of crystallization. However, 
the effective mass of excitonic molecules is 2-3 orders of magni- 
tude smaller than that of the H2 molecule, and crystallization of 
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a biexcitonic liquid probably should not occur at all. Therefore, 
at sufficiently low temperatures a biexcitonic liquid can go into 
a superfluid state. However, as it turns out, the formation of 
a biexcitonic liquid itself is extremely problematical. Keldysh[§| 
and later Brinkman and RicejO] pointed out two important differ- 
ences between H2 molecules and EMs in typical semiconductors. 
First, since the effective masses of the electron and hole are often 
of the same order, the binding energy of the excitonic molecule, 
measured in units of the binding energy of the exciton (atom) , 
should be much less than 0.35 for the hydrogen molecule. Second, 
the contribution of the energy of zero-point motion is consider- 
ably larger in a system of biexcitons because of the small mass of 
a biexciton. For these two reasons, as experiments have shown, in 
typical semiconductors (Si, Ge, and a number of others) at high 
excitation intensities an electron-hole Fermi liquid is formed. One 
therefore expects that the formation of a biexcitonic Bose liquid 
can occur only in crystals in which the effective masses of the elec- 
tron and hole are sufficiently different. Furthermore, the mass of 
the EM should be sufficiently large that the contribution of the 
zero-point motion is sufficiently small and it becomes possible for 
a molecular liquid to form. At the same time, the mass of the EM 
must be small enough that "early" crystallization cannot occur. 

Crystals of monoclinic zinc diphosphide (/5-ZnP2) meet the 
stated requirements. This crystal is characterized by an appre- 
ciable difference of the electron and hole effective masses {a = 
Tne/nih = 0.06), a rather large value of the translational mass 
of the biexciton mbex = 3mo (mo is the free electron mass), and 
a very strong anisotropy of the effective masses of the carriers. 
For example, the anisotropy parameter 7 = f^w/fi^ is equal to 

0. 19 (;U-|| and fi± are the values of the reduced masses of an ex- 
citon in the directions parallel to and perpendicular to the Z 
(c) axis of the crystal). As we know, the anisotropy leads to an 
increase in the binding energy of free electron-hole complexes, 

1. e., excitons, biexcitons, and their condensates. Therefore, the 
excitonic molecule in /3-ZnP2 has a rather high binding energy 
El^ = 6.7 meV = O.lbE^. 

The excitonic state with the lowest energy in this crystal is the 
dipole-forbidden IS state of the orthoexciton. [jT^ This makes it 
rather easy to create an appreciable concentration of excitons and. 
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hence, biexcitons under laser excitation. Studies by the authors 
have shown that the condensation of biexcitons in /3-ZnP2 crystals 
occurs via a hydrogenlike scenario, i.e., unlike the case of typical 
semiconductors, a liquid of the molecular (insulator) type forms in 
them. In the photoluminescence (PL) spectrum of these crystals 
one observes the so-called C line, which is due to two-photon 
annihilation of biexcitionsQ condensed into a molecular (insulator) 
hquid.[|ri], [12], |T3| In Refs.[|T2] and [tT3j the phase diagram of a 
biexcitonic liquid was measured and its critical parameters were 
determined: Tc = 4.9 K, nbex,c = 4.1 x 10^^ cm"^ (?^bex,c ~ 
63 A = 4.2 aex, where Oex is the excitonic Bohr radius). 

In this paper we propose a model for the quantum biexcitonic 
Bose liquid and calculate its emission spectrum. In the framework 
of this model we analyze the experimentally observed features of 
the fine structure of the emission line of a biexcitonic liquid in 
/3-ZnP2 (the C line) as the temperature and the intensity of the 
laser excitation of the crystal are varied. 



Dispersion relation for elementary excitations of a quan- 
tum biexcitonic liquid 

Thus we assume that the biexcitonic liquid does not crystallize 
down to the temperatures at which quantum effects become im- 
portant. What sorts of elementary excitations can exist in a quan- 
tum biexcitonic Bose liquid? It is logical to assume that, first, 
there are acoustic phonons. Their dispersion relation is given as 

E{k) = nuk, (1) 

where u is the sound velocity, and k is the wave vector. Second, we 
assume that owing to the intermolecular interaction in the liquid 
a collectivization of the intramolecular excitations can occur in 
it. Since the lowest-energy excited state of an excitonic molecule 
is a rotational state, [|T5j these collectivized molecular excitations 



can be rotational excitations of the molecules. Here, following 

^ At two-photon annihilation there is no recoil particle that could take an appreciable part 
of the momentum of the annihilating EM. For this reason, since the photon momentum is 
very small, momentum conservation allows two-photon decay only for an EM with fc « 0, i.e., 
the two-photon PL line should be very narrow, jl^ In two-electron (one-photon) transitions, 
in which only one exciton of a molecule annihilates, the second, surviving exciton can take 
an arbitrary quasimomentum, i.e., the one-photon PL lines are rather wide. This makes it 
possible to distinguish the two-photon PL lines of free biexcitons and of a biexcitonic liquid, 
unlike the case of the one-photon PL lines. 



4 



Landau, [m these collectivized excitations/quasiparticles will be 
called rotons. Rotons in a biexcitonic (molecular) liquid can be 
regarded as molecular Frenkel excitons. The excitation energy 
of the molecule in the presence of phase correlation between the 
molecules in the liquid can be written as 



AE = AEq + D + 2M cos[(A; - k'^ 



)a 



where AE'o is the excitation energy of a free molecule, D is the 
change of the interaction energy of a given molecule with its neigh- 
bors under its excitation, M is the matrix element for the transfer 
of the energy of excitation from the excited molecule to a neigh- 
bor being in the ground state, a is the average distance between 
these molecules, and k is the wave vector of a Frenkel exciton, 
i.e., in our case the roton. The last term in Eq. (2) describes 
the dispersion relation of rotons and is physically meaningful for 
I A A; I = |/c — /cqI < 7r/2a. Thus the photon-roton dispersion re- 
lation for elementary excitations of a biexcitonic liquid can be 
written as 

E{k) = hukatO<k<kQ- 7r/2a (3) 
and k > /cq + 7r/2a ; 



E{k) = huk — 7 cos[(/c — /cq 



a 



at /cq — 7r/2a < A; < /cq + 7r/2a , 

where 7 = —2M. For small |AA;| the second relation of system 
(3) can be written in the form 

ZTTlr 

Equation (4) is exactly the same as the equation proposed by 
Landau for describing the dispersion relation of rotons in su- 
perfluid He IL Thus our choice of the term "roton" is not arbitrary 
but is based on the similarity of the dispersion relations of rotons 
in helium and the collectivized molecular excitations of a biexci- 
tonic liquid. In Eq. (4) we have introduced the following notation: 
TJir = —h^/2Ma^ is the roton effective mass, k^ = k^ — rriru/fi is 
the wave vector corresponding to the roton minimum of the dis- 
persion relation, and A = UukQ + 2M + mrV? /2 is the energy 
gap in the spectrum of elementary excitations of the biexcitonic 
liquid. 
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Starting from the dispersion relation (3) for a quantum biex- 
citonic liquid, let us calculate the PL spectrum for such a liquid. 
The spectrum we are looking for can be obtained by convolution 
of the slit function and the function I{E) = lQp{E)f{E), where 
Iq = const, p(E) = pok'^{E){dk{E) /dE) is the density of states of 
the liquid, which is determined from the dispersion relation for the 
elementary excitations (po = const); f{E) = l/[exp(^//csT) — 1] 
is the distribution function for the excitations/quasiparticles of 
the liquid (the Bose-Einstein distribution function). The zero 
of energy E is taken as the value corresponding to the radiative 
transition from the state with k = 0. For the parameters of the 
dispersion curve we used the values obtained from a best fit of the 
calculated luminescence spectrum to the experimental PL spec- 
trum of a biexcitonic liquid, recorded for a /?-ZnP2 crystal of high 
purity. Figures 1, 2, 3, and 4 show the calculated PL spectra of a 
quantum biexcitonic liquid of various densities at different tem- 
peratures and also the characteristics of these spectra as functions 
of temperature T and the square of the excitation intensity I^^^ 
(the density of the biexcitonic liquid is proportional to /exc)- 
is seen that the PL spectrum of a biexcitonic liquid should have 
a two-component structure. In the proposed model it is assumed 
that the shape and parameters of the dispersion curve do not de- 
pend on temperature and depend only on the density of the liquid 
. The density of the hquid influences the shape of the dispersion 
curve and, through it, the shape of the PL spectrum. We used 
the following proportionalities relating the parameters of the dis- 
persion curve and the density of the liquid (which is proportional 
to the square of the excitation intensity): average intermolecu- 
lar distance a oc n~^/^ oc I^^^^^, sound velocity u oc n^/^ oc /exc, 
k'Qtt = const, and the modulus of the matrix element for the exci- 
tation transfer between molecules of the liquid \M\ oc oc I^^^. 

Concluding this Section we must mention the following. A 
short-range order should be established in the system of biexci- 
tons due to the substantial interaction between them; however, 
there is apparently no long-range order. The presence of short- 
range order ensures that relation (2) will apply at least in a qual- 
itative way. 
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Experimental details 



We studied semiconducting single crystals of monoclinic zinc diphos- 
phide of high purity, grown from the gas phase. The technology of 
the synthesis and growth of these crystals is described in Ref. []T^ . 



For excitation of the luminescence we used a cw Ar+ laser (emis- 
sion line 5145 A). 

During laser excitation the temperature of the crystal at the 
point of excitation differs substantially from the temperature of 
the external environment, i.e., of the helium bath surrounding 
the sample and the unilluminated region of the sample. This 
makes it impossible to use a thermosensor for precision tempera- 
ture measurements, as one would get some averaged temperature. 
It was therefore necessary to use an internal temperature standard 
(ITS) that would make it possible to determine the temperature 
of the crystal at the point of excitation. For this standard we 
used the temperature dependence of the spectral position of the 
narrow, intense B line of the PL spectrum (the emission line of 
the free forbidden orthoexciton). Because of the variation of the 
energy gap Eg with variation of temperature, the B line changes 
its spectral position (shifts to lower energy with increasing tem- 
perature). Since the shift of the B line with temperature is small 
[dT/dXB = 3.625 K/A), for correct measurement we recorded a 
reference line of neon in the investigated spectral region simulta- 
neously with the PL spectrum. This technique made it possible 
to determine the temperature of the sample at the point of exci- 
tation to an accuracy of 0.05 K. 



Biexcitonic liquid in /3-ZnP2 crystals. Experimental re- 
sults and discussion 

At excitation intensities above 1 kW/cm^ and temperatures below 
5 K the PL spectrum of /3-ZnP2 crystals, as we have said, contains 
an emission line of the biexcitonic liquid (the C line, henceforth 
called the C spectrum). As we see from Figs. 5, 6, 7, 8, and 9, the 
C spectrum has a two-component form, as was predicted by the 
model set forth above. Using the dispersion relation of elementary 
excitations of a biexcitonic liquid (3) we fit the experimental PL 
spectrum of the biexcitonic liquid (Fig. 5), which was obtained for 
a sample of high purity (excitation energy 3 kW/cm^). This made 
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it possible to determine the parameters of the dispersion curve of 
the elementary excitations: A = 0.49 meV, ko = 9.8 x 10^ cm~-^, 
w = 1.4 X 10^ cm/s, and rrir = 2.2mo (mg is the free electron 
mass), and also the value of the intermolecular distance in the 
liquid a = 92 A = G.lagx and the temperature of the crystal at the 
point of excitation, T = 1.5 K. The values obtained are physically 
reasonable, a fact which, in our view, tends to confirm that the 
proposed theoretical model is correct. Further evidence of this is 
the good agreement of the experimental and theoretical spectra, 
in view of the approximate, model character of the dispersion 
curve given by Eqs. (3). 

The evolution of the shape of the C spectrum and the behav- 
ior of its characteristics under variation of temperature and of the 
square of the excitation intensity are shown in Figs. 6, 7, 8, and 
9. The behavior of the experimental spectrum under variation of 
temperature (Figs. 6 and 7) is similar to that of the theoretical 
spectrum (Figs. 1 and 2). Analysis of the C spectrum under vari- 
ation of the excitation intensity is made difficult by the fact that 
simultaneously with the variation of excitation intensity there is 
also a variation of temperature, which has a definite effect on the 
shape of the spectrum. In comparing Figs. 2, 4, and 9 we can 
conclude that the change in the C spectrum at moderate levels of 
excitation is probably due mainly to the influence of temperature 
(the rise in temperature with increasing /exc)? while at higher ex- 
citation intensities the evolution of the shape of the C spectrum 
occurs mainly on account of the increase in the density of the 
liquid with increasing /exc- What we have said agrees with the in- 
tensity dependence of the temperature of the crystal at the point 
of excitation, shown in Fig. 10. Thus the experiment is described 
by the theory in a completely satisfactory way. The shape of the 
C spectrum can differ substantially for different samples. This is 
due to the fact that at the same values of the excitation intensity 
and crystal temperature, a liquid of higher density should arise in 
purer samples than in less pure samples. Consequently, the shape 
of the PL spectrum of the liquid varies with its density. 

The results presented above suggest that a quantum biexci- 
tonic liquid characterized by a phonon-roton dispersion relation 
for the elementary excitations can form in /?-ZnP2 crystals under 
certain conditions. This sort of dispersion relation, as we know. 
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is evidence that superfluidity can arise in a Bose liquid, in par- 
ticular, in liquid ^He. Is it possible that superfluidity can arise in 
a biexcitonic liquid in a /?-ZnP2 crystal? To answer this question 
we performed the following experiments. For some of the sam- 
ples (Fig. 10) we studied the temperature of the crystal at the 
point of excitation as a function of the intensity of the excitation. 
The normal situation is for the temperature of the crystal to in- 
creasing with increasing excitation intensity. Besides the normal 
monotonically increases behavior of T(I^^) (curve 1), for several 
samples we obtained anomalous T{I^^) curves: decreasing (curve 
2) or nonmonotonic (curve 3). We think that the anomalous 
T[ll^ dependences can be explained by an anomalously large 
(giant) thermal conductivity, which accompanies the appearance 
of the superfluid (s) component in a quantum liquid below the 
superfluid transition. Such an effect is well known for superfluid 
He II. We propose the following explanation for the anomalous 
behaviors observed. At a certain excitation intensity the density 
of the liquid reaches values sufficient for a transition of the liquid 
to the superfluid state. Upon further increase in /exc the density 
of the liquid increases, with the result that the temperature of 
the superfluid transition increases and, with it, the fraction of 
the s component. This can increase the thermal conductivity, 
i.e., the thermal energy should be removed more efficiently from 
the excitation region of the crystal, and so the temperature of 
the crystal should decrease with increasing excitation intensity. 
Apparently those crystals with a monotonically decreasing T(I^^) 
curve (curve 2 in Fig. 10) are the purest and most perfect: they 
have a low concentration of lattice defects, and therefore a biex- 
citonic liquid is formed in them with a sufficiently high density 
for the onset of superfluidity and the corresponding giant thermal 
conductivity. Those crystals with a nonmonotonic T{I'^^) curve 
(curve 3 in Fig. 10) would be less pure. The liquid formed in 
them is of a lower density but is nevertheless sufficient for the 
appearance of the s component. However this component repre- 
sents a smaller fraction than in the perfect crystals, and the giant 
thermal conductivity provides a sufficiently effective heat removal 
only up to certain values of the incoming laser power. When the 
power is increased further, the temperature of the crystal begins 
to grow. In those crystals with a monotonically increasing T(/g^) 
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dependence (curve 1 in Fig. 10) the concentration of defects is 
rather high, and therefore the density of the hquid in them does 
not reach the values necessary for a transition to the superfluid 
state and the resulting onset of giant thermal conductivity. As a 
result, the temperature of the crystal at the point of excitation 
increases with increasing excitation intensity. 

This study was supported in part by the Foundation for Basic 
Research of the Ministry of Scientific and Technological Affairs of 
Ukraine (Project No. 2.4/311). 
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Figure captions 



Figurel. Evolution with temperature of the calculated pho- 
toluminescence spectrum of a quantum biexcitonic liquid with a 
phonon-roton dispersion relation for the elementary excitations. 
The arrow labeled / indicates the low-energy (/) component of the 
spectrum, and the arrow h the high-energy (h) component. 

Figure2. Temperature dependence of the characteristics of 
the theoretical PL spectrum of a quantum biexcitonic liquid: 
Ii I /tot is the ratio of the intensity of the I component to the total 
intensity of the C spectrum, Ih/hot is the same for the h com- 
ponent (a); Ii/Ih is the ratio of the intensities of the / and h 
components (b); the positions of the maxima of the components 
(c); the spectral distance (splitting) between the maxima of the 
components (d); the half- widths of the components (e) and (f). 

Figures. Evolution of the theoretical PL spectrum of a quan- 
tum biexcitonic liquid under variation of the square of the exci- 
tation intensity (density of the liquid): n oc ll^. The curves are 
labeled with the value of /^^ in arbitrary units. 

Figure4. Characteristics of the theoretical PL spectrum of 
a quantum biexcitonic liquid as functions of the square of the 
excitation intensity (density of the liquid). The notation is the 
same as in Fig. 2. 

Figures. Experimental PL spectrum of a biexcitonic liquid 
in a /3-ZnP2 crystal (solid curve) and its approximation by the 
theoretical spectrum (dashed curve) (a); the dispersion curve for 
the elementary excitations of the quantum biexcitonic liquid with 
the parameters obtained by fitting the theoretical PL spectrum to 
the experimental spectrum (b). Excitation intensity 3 kW/cm^. 

FigureG. Evolution of the experimental PL spectrum of a 
biexcitonic liquid in a /3-ZnP2 crystal under variation of tempera- 
ture of the sample at the point of excitation. Excitation intensity 
4.8 kW/cm^. 

FigureT. Temperature dependence of the characteristics of 
the experimental PL spectrum of a biexcitonic liquid in a (3- 
ZnP2 crystal. The notation is the same as in Figs. 2 and 4. 
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Plotted along the vertical axis in part Fig. 7c is the difference 
[E^{T) - Eb{T)] - [Em{To) - Ebin)], where Em is the position 
of the maximum of the I or h component, and E^ is the position 
of the maximum of the narrow line of the forbidden free orthoex- 
citon {B line), which was used to determine the temperature of 
the crystal at the point of excitation, and Tq is the minimum 
temperature reached in the experiment. 

Figures. Evolution of the experimental PL spectrum of a 
biexcitonic liquid in a /?-ZnP2 crystal under variation of the square 
of the excitation intensity. The curves are labeled by the value of 
II, in (kW/cm2)2. 

Figure9. Characteristics of the experimental PL spectrum of 
a biexcitonic liquid in a /?-ZnP2 crystal under variation of the 
square of the excitation intensity. The notation is the same as in 
Figs. 2 and 4 and 7. Plotted along the vertical axis in Fig. 9c is 
the difference [^^(4) - ^6(4)] " [EUli) " where E^ 

is the position of the maximum of the I or h component, is 
the position of the maximum of the B line, and /q is the lowest 
excitation intensity used in the experiment. 

FigurelO. Temperature of the crystal versus the square of 
the excitation energy for different samples: ZJ09-9 (i), ZJ09-14 
(^), and ZD99-1 (5). 
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